Oscillating Sign of Drag in High Landau Levels 
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Motivated by experiments, we study the sign of the Coulomb drag voltage in a double layer system 
in a strong magnetic field. We show that the commonly used Fermi Golden Rule approach implicitly 
assumes a linear dependence of intra-layer conductivity on density, and is thus inadequate in strong 
magnetic fields. Going beyond this approach, we show that the drag voltage commonly changes sign 
with density difference between the layers. We predict that in the Quantum Hall regime the Hall 
and longitudinal drag resistivities are comparable. Our results are also relevant for pumping and 
acoustoelectric experiments. 
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Drag experiments in coupled two-dimensional electron 
systems provide information on the response of a system 
at finite frequency and wavevector and are thus comple- 
mentary to standard DC transport measurements KB- 
In a typical drag experiment, a current is applied to the 
active layer of a double-layer system and the voltage Vd 
induced in the other passive layer is measured, with no 
current allowed to flow in that layer. In a simple picture 
of drag, the current in the active layer leads - via inter- 
layer Coulomb or phonon interaction - to a net transfer 
of momentum to the carriers in the passive layer. At 
conditions of zero current in the passive layer, this mo- 
mentum transfer is counteracted by the build-up of the 
drag voltage Vd- In the cases of two electron layers or two 
hole layers, the drag voltage points opposite to the volt- 
age drop in the active layer. This is defined as positive 
drag. Negative drag occurs in systems with one electron 
layer and one hole layer [^]|| . 

Prior theoretical work on drag J|,HJ| was often based 
on or reduced to a Fermi Golden Rule analysis (see, e.g., 
Zheng and MacDonald fl). In this analysis, the sign of 
drag does not vary with magnetic field B, temperature 
T, or difference in Landau level (LL) filling factor v be- 
tween the two layers. By contrast, recent experiments at 
large perpendicular B [|]-|| observe that the sign of drag 
changes with all of these parameters in systems of two 
coupled electron layers - specifically in the Shubnikov- 
de-Haas (SdH) and integer quantum Hall (IQH) regimes. 
Feng et al. || find positive drag whenever the topmost 
partially filled LLs in both layers are either less than half 
filled or more than half filled. Negative drag is observed 
when the topmost LL is less than half-filled in one layer 
and more than half-filled in the other. Even more sur- 
prisingly, Lok et al. Q| maintained that the sign of drag 
was sensitive to the relative orientation of the majority 
spins of the two layers at the Fermi energy. 

In this paper, we show that the use of the Fermi Golden 
Rule approach is inappropriate in the SdH and IQH 



regimes and that a more careful analysis opens different 
routes to drag with changing sign. Remarkably, we find 
that Hall drag can be of the same magnitude as longitu- 
dinal drag in these regimes. We emphasize that a naive 
rationale for the experimental results of Feng et al. 
which points to the similarity between a less (more) than 
half-filled Landau level and an electron (hole) like band 
- leaves out essential physics of the problem. Finally, we 
compare our results to the existing measurements ||[7]] 
and propose some interesting experiments. 

We start from a general linear-response expression 0] 
which relates the drag conductivity a D to the rectifica- 
tion coefficients r(q, ui) of the active and passive layers, 
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The drag conductivity gives the current response J p in 
the passive (dragged) layer via Jf — cfjEj to an applied 
electric field E a in the active layer. Its relation to the 
commonly measured drag resistivity is clarified below. 
In Eq. (Q), U(q,ui) is the screened interlayer (Coulomb 
or phonon-mediated) interaction and the rectification co- 
efficient r(q, u) is defined by 
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In Eq. (||), 3 dc is the DC current induced in quadratic 
response to the driving force exerted by a screened po- 
tential 0(q, u>) of wave-vector q and frequency u>. 

The physical interpretation of the drag expression Eq. 
(|l|) is that the voltage in the active layer creates an asym- 
metry (between q and — q) in the thermal density fluc- 
tuations in that layer. These fluctuations are transferred 
to the passive layer - via Coulomb or phonon interaction 



1 



— where they are rectified to create a current. In experi- 
ment, usually the drag resisitivity p D is measured, which 
gives the electric field response E p in the passive layer to 



a current 3 a driven in the active layer via Ef = pfjJf ■ 
For weak drag, it is related to the drag conductivity by 



P D = p?a D p° 
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where pPW is the resistivity tensor of the passive (active) 
layer. 

A simple but instructive approach to rectification con- 
siders situations in which current and field are locally 
related, J(r,t) = — a(n(t,t))V^>(r,t), and the conductiv- 
ity depends on position and time only through the lo- 
cal density n(r,t) ||. In addition to the current, the 
perturbation </>(q, uj) also induces a density perturba- 
tion (5n(q, uj) = II(q, w)e0(q, w) due to the polarizability 
Il(q, uj). Up to quadratic order in the applied potential, 
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Taking the time and space average of the second term 
then yields a DC contribution to the current given by 

Jdc = Y^ (^) [ n (q^) e< Mq: w )"#(-q> -<*>)] > ( 5 ) 



or equivalently the rectification coefficient is given by 
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At zero magnetic field, the conductivity is local when 
considering scales large compared to the mean free path 
£ c \, implying that Eq. wi) holds in the diffusive regime 
defined by ql e \ <C 1 |[I0| . At finite B, and for short 
range disorder, Eq. (g) is still valid in the diffusive regime 
WT,Dq 2 T <C 1, where D and r denote the appropriate 
diffusion constant and scattering time. Indeed, we have 
checked this ]ll[ | by an explicit diagrammatic calculation 
in the self-consistent Born approximation (SCBA) Jl2| . 
Specifically, once the magnetic field is strong enough such 
that the cyclotron radius R c = Tik F /(eB) is small com- 
pared to £ e l and for short-range correlated disorder, the 
conductivity is local on scales larger than R c , and Eq. 
(g) holds for qR c < 1 |[(J. Since R c <C 4i even for very 
small magnetic fields, and since the diffusion constant de- 
creases rapidly with the application of a magnetic field, 
Eq. (g) is expected to be applicable for a significant part 
of the wave- vector range of the integration in (|l|) for typ- 
ical experiments. 

Using Eq. in the general drag expressions, Eqs. (|l|) 
and (g), one readily finds that up to an overall positive 
prefactor, the drag resistivity tensor becomes 



, da p da a 
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This expression easily reproduces some standard results. 
In the absence of a magnetic field, the tensor struc- 
ture is trivial. Observing that the conductivity increases 
(decreases) with increasing electron density for electron 
(hole) layers, we recover that coupled layers with the 
same type of charge carriers exhibit positive drag, while 
coupled electron-hole layers have a negative drag resis- 
tivity. If the system is strictly electron-hole symmetric, 
the derivative of the conductivity with respect to density 
vanishes and consequently, there is no drag. Treating the 
conductivity in the presence of a magnetic field in a sim- 
ple Drude-type picture, a depends linearly on the density 
en (so long as the scattering time is taken to be density 
independent). Eq. (g) then reproduces results of previous 
works [QMg] t na t calculate drag in a Fermi Golden Rule 
or scattering time approximation. In particular, in this 
approximation it is immediately found from Eq. (g) that 
the drag resistivity is diagonal and Hall drag vanishes 
identically. Finite Hall drag can appear even in a Drude- 
type approximation once the scattering time is taken to 
be energy and thus density dependent ff3|| . 

Interesting new effects appear in the presence of a 
strong magnetic field, in the SdH and IQH regimes, where 
the derivative of o~ xx changes sign as the magnetic field or 
Fermi energy is varied. While this is superficially quite 
reminiscent of the experimental results of Refs. [g,f7|, the 
details can be involved due to the many terms that con- 
tribute, once the tensor products in Eq. (g) are multiplied 
out. In the experimental samples, typically p xy 3> p xx 
already at very small magnetic fields. Thus, up to an 
overall positive prefactor 
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Generally the derivative of the longitudinal conductiv- 
ity o xx changes sign in the SdH and IQH regimes, being 
positive for less than half filling of the topmost LL and 
negative for more than half filling. By contrast, the Hall 
conductivity generally increases monotonically with n, 
and its derivative is therefore positive. 

It is not clear, apriori, which of the two terms in 
the curly brackets of Eq. (||) dominates. We obtain 
an oscillatory sign of drag, similar to the experimen- 
tal results in Ref. H when do~ xx /d(eri) dominates over 
do~ xy /d{en). However, if that were the case, drag would 
be negative for equal filling of the two layers (because 
p xy = —p y x)i in contrast to the experimental observa- 
tions. Positive drag for equal filling is only obtained 
when do~ xy /d(en) 3> do~ xx /d(en). Were that the case, 
however, there would presumably be no sign changes of 
the drag resistivity. In fact, in the IQH regime the deriva- 
tives of both components of the conductivity tensor are 
experimentally of the same order (both change by ap- 
proximately e 2 /h in the region of the plateau transition). 
While neither of the two limits is therefore realized in 
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experiment, it is clear that Eq. (@) makes it difficult to 
obtain both positive drag for equal densities and a drag 
sign that oscillates with the difference in densities. 

A striking consequence of Eq. (Jfy is that if we assume 
that da xx j d(en) and da xy / d(en) are of comparable mag- 
nitude in the integer-quantum-Hall regime then the Hall 
drag is of the same order as longitudinal drag. In fact, 
we have 
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Moreover, Hall drag generally changes sign with filling 
factor difference between the two layers. 

Although the experimental disorder potential is long 
ranged, it is instructive to consider also the case of 
short-range disorder, for which the SCBA |l2| becomes 
exact in the limit of high LLs H]. Assuming that 
the disorder broadening is small compared to the LL 
spacing, the conductivity tensor a becomes Q a xx — 
{e 2 /ir 2 h)N{l - [{E F - Sat)/2 7 ] 2 } and a xy = ~{en/B) + 
(e 2 /7T 2 h)(2 1 /nu Jc )N{l - \{E F - E N )/2 1 } 2 } 3 / 2 where 7 
denotes the LL broadening, and the Fermi energy Ef 
is assumed to lie within the iV-th LL (of energy Em = 
(N + \/2)hu> c in the clean case). Thus \da xx /d(en)\ 3> 
\da xy /d(en)\ and da xx /d(en) changes sign as a function 
of filling - being positive for less than half filled LLs and 
negative for more than half-filled LLs. Neglecting the 
derivatives of the Hall conductivity, one obtains for the 
diagonal drag resistivity of an isotropic system 
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with C xx a positive function. Thus, we again find the 
surprising result that drag is negative for equal densities 
in the two layers (since \p xy / p xx \ > tt)- More generally, 
drag oscillates with difference in density between the lay- 
ers. It is negative whenever the topmost occupied LLs in 
the two layers are both less than half filled or both more 
than half-filled. It is positive if the topmost LL is less 
than half-filled in one layer and more than half-filled in 
the other. Moreover, the Hall drag resistivity is compa- 
rable to the diagonal drag resistivity also in this model 
situation. In fact, one has 



D 
Pxy 



d(en) d(en) 



r p a I 
[rxyrxx * 



jp a "i 
rxxPxy] 



(11) 



with C xy = C. 



Unlike before, this result is now a 
consequence of the fact that p xx and p xy are of similar 
magnitudes. 

When the LL broadening is comparable to but exceeds 
the LL spacing (lu c t < 1), we find [ [TT| from a numerical 
evaluation of the SCBA equations that by contrast, the 
derivative of a xy dominates over that of a xx while still 
changing sign as function of filling (for Nj/hu> c 3> 1). 
However, in this limit, a xx dominates over a xy (and hence 



Pxx over Pxy) so that again, we find negative drag for 
equal filling of both layers. 

Experiments do not necessarily satisfy the condition 
Dq 2 r <C 1 for the diffusive regime. Thus, we now turn 
to a discussion of drag in the ballistic regime where 
Dq 2 r » 1. At zero B, Eq. (|) holds also in the bal- 
listic regime @. Indeed, Eq. (§) with <r given by the 
Drude expression is also derivable within a Boltzmann 
approach which should be valid for q/k F ,u>/ E F -C 1 and 
at fields low enough that SdH oscillations are absent . 
Surprisingly, we find that this is not generally true for 
higher magnetic fields. While we have not succeeded in 
deriving a general expression, analogous to Eq. (|]), in 
this regime, we computed the rectification explicitly in 
the high-magnetic field limit introduced above. 

Our calculation starts from the diagrammatic expres- 
sion for the rectification coefficient S 

r(q,w) = — Tr{lG-e iqr [G- - G+]e- lqr G + ) (12) 

7TZ 

Here, G^ denotes the impurity-averaged Green function 
in the SCBA, I is the current operator, and Tr the trace 
over the single-particle states. The calculation is simpli- 
fied by the fact that vertex corrections can be neglected 
in the ballistic limit. 

Using standard results for the matrix elements of the 
current and density operators between LL wavefunctions 
and exploiting the small parameter ~f/fwj c <C 1, we ob- 
tain for the longitudinal rectification 



r,, = 



-AwR t 
N 



Jo(qRc)Jl(qRc) 



dn \ da r 



dp J d(en) 
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with Ji the i th Bessel function and N the LL number. A 
similar calculation for the transverse rectification gives 
L ± = (l/w c r)ri, or 



2 (-AuR, 
r± = 3 



N 



MqR c )J l{ qR c )( d -f] pL. (14) 
\dp J d(en) 

While both T\\ and T± still include factors that can be 
written as derivatives of the corresponding conductiv- 
ity, similar to the diffusive regime, the prefactors are no 
longer equal to one another and cannot be expressed in 
terms of the polarization operator Tl(q,uj). Remarkably, 
the rectification in the ballistic limit also changes sign 
with qR c due to the matrix elements of the density, in 
addition to the sign changes of do j d(en) discussed above. 
Note that such behavior would be impossible in Eq. (|^) 
since the sign of Imn, being fixed by the sign of w, does 
not change with q. 

It is an interesting question whether these additional 
sign changes are reflected in the drag. This is not obvious 
since the oscillating Bessel functions depend on q which 
needs to be integrated over to obtain the drag conductiv- 
ity (cf. Eq. ([!])). For Coulomb drag, the screened inter- 
layer interaction can be approximated by the Thomas- 
Fermi result t/(q, to) — ire 2 q/(K a K p smhqd) with K a (k p ) 
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denoting the Thomas-Fermi momenta of the active (pas- 
sive) layer and d the interlayer distance. Thus, the lo and 
q integrations in Eq. (|l|) factorize. One readily observes 
that the q integration is not dominated by the upper limit 
in the relevant region 1/R C <C q <C l/d, and thus there 
are no additional oscillations from the matrix elements in 
Coulomb drag. This result is also confirmed by straight- 
forward numerical evaluation. 

For phonon drag, the situation is slightly different. 
Using a qualitative approximation ||, \U(q, uj)\ 2 as 
Kq 3 5(u> — qv) with v the sound velocity, the integrals 
in Eq. (0) yield an additional factor T A [g[q T \R a c - BP C \) + 
g{q T \R a c + R p c \)] where q T = {T/T B )k F with T B = 2hvk F 
the Bloch-Griinciscn temperature (for typical samples 
Tb ~ 10 K). Here, g(x) is a smooth function (which can 
be calculated analytically) that starts with g(0) — 1 and 
decreases for small x crossing zero at x ~ .04 then reach- 
ing a minimum of g(.075) w —.37 then increasing mono- 
tonically and exponentially towards zero. For all but the 
lowest temperatures, the term involoving \R^ — R^\ domi- 
nates, and thus the overall sign of the drag changes when 
q T \R a c - BP C \ = \v a - v p \T/T B w .04. 

In passing, we note that while the sign changes of 
r(q, u>) with qR c in the ballistic limit seem to have 
only minor consequences for drag, they should be ob- 
servable in other experiments of recent interest, namely 
in pumping | fl5f and in acoustoelectric experiments 0, 
which also measure the rectification coefficient r(q, lo) 
|jl6|| . In pumping experiments, a potential <fi(r,t) is ap- 
plied to the system which is typically made up of the 
sum of two potentials oscillating out of phase tj>(r, t) = 
4>i(r) cos(uit) + 02(r)cos(wi + 8). Due to the symme- 
tries of r(q, co), the pumping current density is given 
by J dc = (l/2)sin(5)E q r(q,w)Im[^(q)0 2 (-q)]. In 
acoustoelectric experiments, acoustic waves sent through 
a (piezoelectric) crystal apply an "external" electric po- 
tential <f) ext (q, lo) to the electrons in the system with 
lo = cq and c the wave velocity. The driven electric 
current is given precisely by Eq. (|2|) where </>(q, to) = 
4> ext (q,ui)/[l — v(q)U(q,co)] is the screened potential as- 
sociated with <fi ext (v(q) denotes the bare Coulomb inter- 
action). We note that the acoustoelectric experiments of 
Rcf. J^], which were carried out in the diffusive regime, 
have been successfully analysed using Eq. (||). 

For the cases considered here, an oscillating sign of 
drag is always accompanied by negative drag for equal 
fillings. This apparent conflict with the experimental re- 
sults of Feng et al. might be resolved in one of the fol- 
lowing ways: (a) For long-range disorder and hu> c compa- 
rable to the LL broadening, the Hall conductivity might 
become non-monotonous, cf. Eq. ([?]) and the SCBA re- 
sults, (b) In measurements, the large Hall drag resistivity 
might mix into the measured drag response. On the other 
hand, a simple-minded extension of our results to include 
the two spin directions definitely fails to account for the 



apparent spin dependence observed in Ref. 0. We take 
this failure as an indication that the observations of Rcf. 
cannot be explained by theories where correlations 
induced by Coulomb interactions are neglected. 

In conclusion, we have considered drag in high LLs 
with an emphasis on the sign of the effect. Remark- 
ably, we find that both in the ballistic and the diffu- 
sive regimes, drag in high LLs cannot be described by a 
widely-used Fermi Golden Rule expression [Q . Our anal- 
ysis naturally opens the possibility of sign changes of the 
drag resistivity as function of filling factor difference be- 
tween the two layers. Moreover, it implies that Hall drag 
can be of the same order as longitudinal drag in high 
LLs. Surprisingly, we find that the sign of drag can be 
quite different from naive expectations. In particular, in 
several regimes drag is negative for two identical electron 
layers. We believe that it would be particularly inter- 
esting to check experimentally our predictions that Hall 
drag can be of the same order as longitudinal drag and 
that the acoustoelectric current in the ballistic regime 
changes sign as function of the wavevector. 
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